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■^j- . Abstract 

O . 

We consider space and time dependent fuzzy spheres S 2p arising in Dl — D(2p + 1) 
intersections in IIB string theory and collapsing D(2p)-branes in IIA string theory In the 
case of S 2 , where the periodic space and time-dependent solutions can be described by 
Qh! Jacobi elliptic functions, there is a duality of the form r to - which relates the space and 
time dependent solutions. This duality is related to complex multiplication properties of 
the Jacobi elliptic functions. For S 4 funnels, the description of the periodic space and time 
dependent solutions involves the Jacobi Inversion problem on a hyper-elliptic Riemann 
surface of genus 3. Special symmetries of the Riemann surface allow the reduction of the 
problem to one involving a product of genus one surfaces. The symmetries also allow a 
generalisation of the r to - duality. Some of these considerations extend to the case of 
the fuzzy S 6 . 
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1 Introduction 



Fuzzy spheres of two, four, six dimensions arise in a variety of related contexts. On the 
one hand they describe the cross-sections of fuzzy funnels appearing at the intersection of 
Dl-branes with D3, D5 or D7-branes of Type IIB string theory PUI21 El- I n this context it 
is of interest to follow the spatial evolution of the size r of the fuzzy sphere as a function 
of the co-ordinate a along the D-string. At the location of the higher dimensional brane, 
the cross-section of the funnel blows up. These equations for the funnel which arise either 
from the D-string or the D(2p + l)-brane worldvolume, can be generalised to allow for 
time dependence as well as spatial dependence. The purely time-dependent solutions are 
also relevant to the case of spherical bound states of -DO and D2p-branes of Type IIA 
string theory. 

In the case of the fuzzy 2-sphere, there are purely spatial and purely time-dependent 
solutions described in terms of Jacobi elliptic functions. The spatial and time-profiles are 
closely related and the relation follows from an r — > - duality. It is natural to introduce 
a complex variable U\ = a — it. For solutions described in terms of elliptic functions, the 
inversion symmetry is related to the property of complex multiplication ui — > iu\. The 
periodic spatial solutions describe a configuration of alternating branes and anti-branes. 
At the location of the brane or antibrane, the radius r of the funnel blows up. This is 
a well-understood blow-up, expected from the geometry of a 1-brane forming a 3-brane. 
The periodic solutions in time describe collapse followed by expansion of 2-branes. The 
collapse point is a priori a much more mysterious point, where the size of the fuzzy sphere 
is sub-stringy. Nevertheless the r — > - duality following from the equations of the Born- 
Infeld action imply that the zeroes of the time evolution are directly related to the blow-up 
in the spatial profile. 

In the case of the fuzzy 4-sphere, the functions defining the dynamics are naturally 
related to a genus 3 hyper-elliptic curve. Using the conservation laws of the spatial or time 
evolution, the time elapsed or distance along the 1-brane can be expressed in terms of an 
integral of a holomorphic differential on the genus 3 hyper-elliptic curve. The upper limit 
of the integral is the radius r. Inverting the integral to express r in terms of U\ = o — it 
is a problem which can be related to the Jacobi Inversion problem, with a constraint. 
Because of the symmetries of the genus 3 curve, it can be mapped holomorphically to a 
genus 1 and a genus 2 curve. The genus 2 curve can be further mapped to a pair of genus 1 
curves. The Jacobi inversion problem expressed in terms of the genus 2 variables requires 
the introduction of a second complex variable «2 and we find that there is a constraint 
which relates w 2 to U\. As a result, an implicit solution to the constrained Jacobi inversion 
problem can be given in terms of ordinary (genus 1) Jacobi elliptic functions. The solution 
is implicit in the sense that the constraint involved is transcendental and is given in terms 
of elliptic functions. We give several checks of this solution, including a series expansion 
and calculations of the time of collapse or distance to blow-ups. The symmetries which 
allow the reduction of the problem to one involving lower genus Riemann surfaces also 
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provide dualities of the type r — > - which relate poles to zeroes. 

In section 5 we extend some of these discussions to the fuzzy 6-sphere. The space 
and time dependence are related to integrals of a holomorphic differential on a genus 5 
Riemann surface. A simple transformation relates the problem to genus 3. But we have 
not found a further reduction to genus one. The solution r(t, a) can still be related to a 
constrained Jacobi inversion problem, which can be solved in terms of genus 3 Riemann 
theta functions. As far as large-small symmetries are concerned, the story is much the 
same as for S* 4 in the limit of large 'initial' radius tq. In the time-dependent problem by 
'initial' radius we mean the point where the radial velocity is zero. In the spatial problem, 
it is the place where the 4^ = 0. For general r there are still inversion symmetries of the 
type (1 + r 4 ) — > (1 + r 4 ) -1 , but they involve fourth roots when expressed in terms of r, 
so are not as useful. 

Appendix A.l gives a discussion of the BPS condition for the D1-D3 system in a 
Lorentz invariant form appropriate for space and time dependence. We discuss some 
boosts of solutions described in section 2, as well as the relation between the BPS equation, 
the Yang-Mills equation and the DBI equation. Appendix A. 2 uses the Chern-Simons 
terms in the Dl and D3 action, to show that the solutions we have considered do indeed 
carry both Dl and D3 charge. Appendix B describes some Lagrangians related to the 
ones that appear in the Dl — D(2p + 1) system, and which give rise to the equations of 
motion related to higher holomorphic differentials for the Riemann surfaces mentioned 
above. Appendix C describes the derivation of the Jacobi-Cn solution for the fuzzy S 2 
by steps using Weierstrass p functions, since the discussion of section 5 on the fuzzy S 6 
is expressed in terms of higher genus generalisations of the p functions. 

2 Space and Time-Dependent Fuzzy S* 2 

2.1 Non-abelian DBI description of non-static Dl _L D3 funnels 

The static system consisting of a set of N D-strings ending on an orthogonal D3 has been 
thoroughly studied [UI2]- There exist two dual descriptions of the intersection at large- N, 
one from the Dl and one from the D3 worldvolume point of view. In the Dl-picture it 
is described as a funnel of increasing radius as we approach the D3 brane, where the 
D-strings expand into a fuzzy-S* 2 . In the D3-picture the worldvolume solution includes a 
BPS magnetic monopole and the Higgs field is interpreted as a transverse spike. Although 
the Dl picture is valid far from the D3 and the D3 picture close to it, there is a significant 
region of overlap which validates the duality. Here we will enlarge this discussion by lifting 
the static condition. 

We begin by considering the non-abelian DBI action of iV D-strings in a flat back- 
ground and with the gauge fields set to zero 

SdIi = -Tx ! d 2 a STr^-det ( Vab + \ 2 d a &Q~ j 1 d b &) det{Qv) , (2.1) 
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where a, b worldvolume indices, the $'s are worldvolume scalars, A = 2na' = 2ti£ 2 and 

QV = 8*3 (2.2) 
The expansion of this to leading order in A yields the action 

S% l BI ~ —Ti J d 2 a (n + y5Tr ^d a ^d a ¥ + ^[$\ , &]j + (2.3) 

and the following equations of motion at lowest order, which are the Yang-Mills equations 

d a d a & . (2.4) 

We will consider the space-time dependent ansatz 

= R(a, t) a\ i = l,2,3 (2.5) 

where the a*'s are generators of the irreducible N x N matrix representation of the SU(2) 
algebra 

[a\ a j ] = 2ie ijk a k , (2.6) 

with quadratic Casimir Y^=i( al ) 2 = c ^-nxn = (N 2 — l)ljVx7V- The resulting scalar field 
configuration describes a non- commutative fuzzy S 2 with physical radius 

Kh(°> r) = ^Tr[¥(a, r)¥(a, r)\ = ^R 2 . (2.7) 
By replacing the ansatz ()2.5|) into ({23$ we get the non-linear action 



S = - Tl J d 2 a STr\J 1 + X 2 cR' 2 - X 2 cR 2 \J 1 + AX 2 cR 4 . (2.8) 

By varying this with respect to R we recover the full equations of motion. Ignoring 
corrections that come from the application of the symmetrised-trace prescription, which 
are subleading at large N jj], these are given by 

2X 2 cM'k + R"(l - X 2 ck 2 ) - kl + X 2 cR 2 ) = 8R 3 ( 1 + ~ ) . (2.9) 

\ 1 + 4A 2 Cj R 4 J 

We can convert the above formula to dimensionless variables by considering the re- 
scalings 

r=^cR, f=J^f, S=J±*, (2-10) 
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which will then imply 

r < = 4A^, f*y = A^V, f*v =A ^V, (2.ii; 



\df J \ dr I ' \da J \ da 

df* A\dT'') ( ' ' da 2 4\da*) ( ' ' dadf 4 \dadr j 1 ' 
The simplified DBI equations of motion can then be written in a Lorentz-invariant form 

d^r + {d^r) (d u r) {d v r) - (d^r) (d u r) (0"r) = 2 r 3 ( - + ^Iff^ j , (2.12) 

where /i and z/ can take the values <r, f . Further aspects of Lorentz invariance and boosted 
solutions are discussed in Appendix A. 

One can also write down the re-scaled action and energy density of the configura- 
tion, by making use of the dimensionless variables with the dots and primes implying 
differentiation with respect to the re-scaled time and space respectively 



S 3 = - J rfVv/(l + r' 2 -f 2 )(l + r 4 ) 



E = (1 + " 2) 7rS- (2 ' 13) 

We will now switch to the dual picture. The abelian DBI action for a D3-brane with 
a general gauge field, a single transverse scalar in the x 9 direction, in a flat background 
is 

SEbi = ~ T 3 J d 4 a^/-det( Vab + A 2 <9 a $d 6 $ + XF ab ) . (2.14) 
The determinant, for a general gauge field, can be calculated and gives 



Sj>%j = -T 3 J d 4 ay 1 + A 2 |£?| 2 + A 2 (V$) 2 + X 4 (B ■ V$) 2 - A 2 |^| 2 - A 2 $ 2 

-A 4 |£ x V$| 2 - A 4 <j> 2 |5| 2 - A 4 (£ • B) 2 + 2A 4 $ V$ • (B x E) , (2.15) 

from which one can derive the spherically symmetric equations of motion, in the absence 
of electric fields and with a single, radial component for the magnetic field on the D3, 

A$"(l - A 2 $ 2 ) - A$(l + A 2 $' 2 ) + 2A 3 «^ / <i'<i> / = -8^ 3 A$ / (1 . + n A 4 2<l>/2 ~ X *fP . (2.16) 

(4i? D3 + A N ) 

Note that for this configuration, the field $ only depends on the D3 radial co-ordinate 
Rr>3 and time t. Expression ()2.16|) then looks similar to (|2.9jl and we can show that it 
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is indeed the same, when written in terms of Dl world- volume quantities. We consider 
total differentials of the fuzzy sphere physical radius R p h(a,r) 



dR, 



OR 



ph 



ph 



da 



da + 



dR 



ph 



Or 



dr 



(2.17) 



and recover for constant R p h and r respectively 



da 



R 



ph 



R 



ph 



r: 



a$ 



da 



ph 



dR 



ph 



R' 



(2.18) 



ph 



where we are making use of the identifications 1 R p h = Rd3, t = t and a = A$. Then the 
second order derivatives of a are; 



d 



dR, 



ph 



da 
■dr 



A$' 



R 



ph 



d_ 



da 

Rph dT 



A$ 



and 



d 



R 



ph 



dR 



ph 



da 



dR 



ph 



(2.19) 

The Dl-brane solution R p h(a,t) can be inverted to give a(R p h,t). By employing the 
following relations 



df (a(R ph ,t),t) 



dt 



R 



we have 



A$'' 



ph 



JDtt 

^ph 

^ph 



dt 



+ 



df 
da 



da 
t~dl 



and 



df (a(R ph ,t),t) 



R 



ph 



dR 



pit 



and 



A$ 



D/2 

n ph 



RphRph 2RphR P h + Rph j-)/ 
R'L 



= df 


da 


t da 


tdR p h 
( 




TDll \ 

2 ph \ 





(2.20) 



ph 



D/2 
^ph 



R^h Rn 



L ph 1 l ph j^, 



ph 



(2.21) 



ph 



By replacing these into (|2.16|) . one recovers the exact non-linear equations of motion 
(|2.9|) in terms of the physical radius R p h- This guarantees that any space-time dependent 
solutions of (j2.9J) will have a corresponding dual solution on the D3 side. 



2.2 Arrays of branes in space and Collapse/Re-expansion in 
time dependence 

We now restrict to purely time dependent solutions of equation ()2.12|) . The resulting DBI 
equations of motion are identical to those coming from a Lagrangian which describes a 
set of iV DO's, expanded into a fuzzy S 2 . This configuration also has an equivalent dual 
DBI description is in terms of a spherical D2-brane with iV-units of magnetic flux 0j. To 

1 See also Appendix A. 
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simplify the notation, the re-scaled variables r, a of (j2.10j) will be called t, a. Then the 
conserved energy density (j2.13|) (or energy in the D0-D2 context) at large N is 



E = V , . (2.22) 



If ro is the initial radius of the collapsing configuration where r = 0, E = \Jl + r 4 and 

we get 

4 4 

(2.23) 



1 + r, 



This allows us to write 



which can be inverted to give 



dt = I v " , (2.24) 
r Vr 4 - -•' 



r(t) = ±r Cn (t, -^j , (2.25) 

where t = -^tp= . Such solutions were first described in [5] and more recently in jUElHIE] 2 - 

The function r(t) describes a D2-brane of radius starting at r = r , at t = 0. It 
decreases to zero, then goes negative down to a minimum — r and then increases back 
through zero to the initial position. The cycle is then repeated (see Figure 1). The region 
of negative r is somewhat mysterious, but we believe the correct interpretation follows if, 
elaborating on ()2.7jl . we define the physical radius as 

^ = + ^^=W < 2 - 26 > 

in the region of positive r and as 

R rk= "f¥^)^-^ (2-27) 



in the region of negative r. This guarantees that R p h remains positive. The change in 
sign at should not be viewed as a discontinuity that invalidates the use of the derivative 
expansion in the Dirac-Born-Infeld action, since the quantity that appears in the action 
is r (or R) rather than Rph- Continuity of the time derivative <9jr at r = also guarantees 
that the D3-brane (or D2-brane) charge is continuous 3 . This interpretation is compatible 
with the one in [T]. where different signs of the R were interpreted as corresponding to 
either a brane or an anti-brane emerging at the blow-up of the S 2 funnel. 

2 See Appendix C for a derivation of this result using Weierstrass p-functions. 
3 See Appendix A for expressions for the charge. 
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Instead of dropping space dependence we can restrict ourselves to a static problem by 
making our ansatz time independent. There is a conserved pressure T crcr 

By plugging-in the correct expression we recover 



d /1+r 4 , 

, (2.29) 



da V 1 + r' 2 

which can be combined with the initial condition r' = at tq to give 

* = (S) ' (2 ' 30) 

The purely space-dependent and the purely time-dependent equations are related by 
Wick rotation t — > ia. To apply this to the solution ()2.25|) we can use the identity 



c H ix 'VlJ = c^' (2 ' 31) 

which is an example of a complex multiplication formula 9J. Therefore, the first order 
equation for the static configuration has solutions in terms of the Jacobi elliptic functions 

r{a) = ±r -^—-^ , (2.32) 

Cn (a, ^) 

where a = and it can be verified that these also satisfy the full DBI equations of 

V r o +1 

motion. This solution is not BPS and does not satisfy the YM. The relations between the 
DBI, YM and BPS equations are discussed in Appendix A. 

The r(a) plot reveals that it represents an infinite, periodic, alternating brane-anti- 
brane array, with .Dl-funnels extending between them. The values of o where r blows 
up, i.e. the poles of the Cn-function, correspond to locations of D3-branes and anti-D3- 
branes. This follows because the derivative |^ changes sign between successive poles. This 
derivative appears in the computation of the D3-charge from the Chern-Simons terms in 
the Dl-worldvolume. Alternatively, we can pick an oriented set of axes on one brane and 
transport it along the funnel to the neighbouring brane to find that the orientation has 
changed. On the left and right of a blow-up point the sign of the derivative is the same 
which is consistent with the fact that the charge of a brane measured from either the left 
or right should give the same answer. 

This type of solution captures the known results of F and D-strings stretching between 
D3 and anti-D3's [TTT^ fTT] by restricting to a half-period of the elliptic function in the space 



S 



evolution. It is also possible to recover the BPS configurations of [2] which were obtained 
by considering the minimum energy condition of the static funnel, where r = 



d a & = ± % -e i i k [¥, . (2.33) 

This is equivalent to the Nahm equation ^2] and is also the BPS condition. In dimension- 
less variables it translates to r' 2 = r 4 and has a solution in terms of r = ±1 / (a — (Too), with 
(Joe denoting the point in space where the funnel blows-up. We will restrict the general 
solution to a quarter-period and consider the expansion around the first blow-up which 
occurs at Cn(K(^), ^=), i.e. close to the D3-brane. We get 



r 



ro 



'l+rg 



This is of the form r = — — — — which goes to r = — as r — > 0. 



2.3 An r — ► 1/r duality between the Euclidean and Lorentzian 
DBI equations 

Defining s = § the equation (j2~2l?|) can be written as 

4 4 

- hk < 2 ' 35 » 

Viewing s and r as two complex variables constrained by one equation, this defines a 
genus one Riemann surface. The quantity — which gives the infinitesimal time elapsed is 
an interesting geometrical quantity related to the Riemann surface, i.e. the holomorphic 
differential. 

_____ r 2 

The curve (|2.35J) has a number of automorphisms of interest. One checks that R = — , 

2 

5 = -pr leaves the equation of the curve invariant. This automorphism of the curve leads 
directly to the complex multiplication identity (|2.31|) which relates the spatial and time- 
dependent solutions. Another, not unrelated, automorphism acts as R — > i, R — > ^, 
_"=■*§. The relation between s and s is equivalent to a Wick rotation of the time variable. 
The transformation r — ► - can also be taken to act on the second order equation since 
it does not involve ro (which does not appear in the second order equation). The spatial 
BPS solution to the second order equation can be acted upon by this transformation. 
The outcome is a time-dependent solution describing a brane collapsing at the speed of 
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Figure 1: Analytic plot of the Jacobi elliptic function solution for the static fuzzy-5* 2 
funnel array and the collapsing 2-sphere for t-q = 1. 



light r = ±(£ — too). This solution can be derived as a tq — > oo limit of the general 
time-dependent elliptic solution, in much the same way as the BPS solution was derived 
as an r — > limit of the spatial elliptic solutions. 

The action of the r — > - transformation on the second order equations can be seen 
explicitly. In the case of pure time dependence, the equation of motion in dimensionless 
variables is 




(2.36) 



while in the case of pure spatial dependence 




(2.37) 



A substitution r = 1/R can be used to transform ()2.36|) using 



r 



r 
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to get 



R 2 , ^,^2 2 1 — 



-R^m {R) = ^TTW' (2 ' 38) 

which can be simplified to 

So the effect of transforming r — > i and renaming £ — * a is the same as the substitu- 
tion t — > zcr. This explains the relation between ()2.25|) and ()2.32|) which was previously 
obtained using the complex multiplication property (|2.31|) of Jacobi-Cn functions. 



3 Space and Time dependent fuzzy S A 

3.1 Equivalence of the action for the Dl _L D5 intersection 

We will extend the consideration of space and time dependent solutions to DBI for the 
case of the Dl _L D5 intersection, which involves a fuzzy S 4 , generalising the purely 
spatial discussion of The equations are also relevant to the time-dependence of fuzzy 
spherical D0-D4 systems which have been studied in the Yang-Mills limit in |13j . On the 
Dl side, we will have five transverse scalar fields, satisfying the ansatz 

¥{a,r) =±R{a,r)G\ 2 = 1,... 5, (3.1) 

where the G l, s are given by the action of 5*0(5) gamma matrices on the totally symmetric 
n-fold tensor product of the basic spinor, the dimension of which is related to n by 

N = (n + l)(n + 2)(n + 3) (g ^ 

The radial profile and the fuzzy-S* 4 physical radius are again related by 

R p h{cr, t) = y/c\R(a, r) , (3.3) 

where c is the 'Casimir' G l G l = cIn = n(n + 4)1at x at. By plugging the ansatz (|3.1jl 
into the action and by considering the large- N behaviour of the configuration, one gets 
in dimensionless variables defined just as for the S 2 case in ()2.10j) . 



S 1 = -NT t J d 2 aVl + r /2 - r 2 (l + r 4 ) 



(3.4) 



11 



As in the Dl _L D3 case, we can write the equations of motion for this configuration 
in a Lorentz-invariant way. The result is the same as before with the exception of the 
pre-factor on the right-hand-side 

d^r + (d^r) (d v r) (d v r) - (d^r) (d u r) (d"r) = 4 r 3 ^ 1 + j^^^ j , (3.5) 

where again /i and v can take the values a, r. 

Let us now look at the D5 side: The world-volume action for n D5-branes, with one 
transverse scalar excited, is as usual 

S 5 = -T 5 J (faStr^-det(G ab + A 2 <9 a $<9 b $ + \F ab ) . (3.6) 

Introducing spherical co-ordinates with radius Rd5 and angles a 1 (i = 1, ... ,4), we 
will have ds 2 = —dt 2 + dR? D5 + R 2 Db gijdot' ''da?\ where is the metric of a unit four- 
sphere, the volume of which is given by f d 4 a^fg = 8n 2 /3. We take the construction of 
|2j with homogeneous instantons and, keeping the same gauge field, we generalise $(i?D5) 
to $(i?D5,t). We will not review this here but just state that by using the above, we can 
reduce the previous expression to 



S 5 = -T s J d 6 a^gVl + A 2 $' 2 - A 2 $ 2 (nR% 5 + ^iVA 2 ^ , (3.7) 
Then by implementing once again the relations (|2.18jl we recover 



S 5 = ~T 5 Jd 4 aAdtAdR ph ^Jl + ^-^(^nR 4 ph + ^NX 2 ) . (3. 



ph ph 

This can be easily manipulated to yield the following result 

S 5 = —TiN J dr A da^/l + R' 2 h -R 2 ph (l + ^K^j ■ (3.9) 

Using |^ ~ f , which holds for large- iV and by once again employing dimensionless vari- 
ables, this becomes 

S 5 = —TiN J dr A daVl + r' 2 - f 2 (l + r 4 ) , (3.10) 

which agrees with ()3.4j) and can be further simplified to the re-scaled action 

S 5 = - J d 2 aVl + r' 2 - f 2 (l + r 4 ) . (3.11) 

Thus, every space-time dependent solution described by the DBI equations of motion 
from the Dl worldvolume will have an equivalent description on the D5 worldvolume. 



12 



3.2 Solutions for the space-dependent fuzzy spheres: Funnels 

Assuming r — 0, the action ()3.11|) becomes 



£ 5 = - / d 2 aVl+r' 2 (1 + 



r 4 ) . (3.12) 



There exists a conserved pressure T CTfT = P given by 



1 + r 4 

P = ^=^, (3.13) 
VI + r' 2 



which can be solved to give r' 



ih-Y \ (1 -! /•'}- 



cW P 2 P 2 

(2r 4 + r 8 ) - (2r 4 + r\) 



Defining 4^ = 5 this can be expressed as 



^ = (r -rMr -rJ 
[1 + rg^ 



The roots in the equation above correspond to (1 + r 4 ) = (1 + r 4 ) and (1 + r 4 ) = — (1+r 4 ). 
The second possibility gives r\ = —2 — r 4 . By differentiating the pressure, a second order 
differential equation can be derived 

d 2 r Ar 3 (l + r' 2 ) 



da 2 (1 + r 4 ) 

4r 3 (l + r A ) 



(3.15) 



P 2 

An integral formula can be written for the distance along the Pl-brane using ()3.14j) 

da= f dr — r ^ — , (3.16) 

J ro v /(2r 4 + r 8) - (2r 4 + r 8 ) 

where we have taken the zero of a to be at the place r = r where r' = 0. We deduce 
from the integral that there is a finite value of a, denoted as E, where r has increased to 
infinity. This is similar to part of the spatial solution of the D1-D3 system described by 
the Cn function (see Figure 1). We can show that the full periodic structure analogous 
to that of the spatial D1-D3 system follows in the D1-D5 system, by using symmetries of 
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the equations and the requirement that the derivative ^ is continuous. Note that ()3.14D 
is symmetric under the operations 



I r : r{a) 
T s : r(a) - 



-> -r(o-) 
r(a-E) . 



(3.17) 



The branch with r increasing from r to oo, as a changes from to E, can be acted on by 
T2sI a Ir to yield a branch where r increases from — oo to — r over E < er < 2E. Acting 
with T2T,I r gives a branch where r decreases from — ro to — oo over 2E < o~ < 3E. Finally 
a transformation of the original branch by T^I a gives r decreasing from oo to ro over 
3S < a < 4E. These four branches patch together without discontinuity in r or ^ and 
can further be translated by integer multiples of 4E, to give a picture qualitatively similar 
to the D1-D3 case, but now describing Dl funnels between alternating D5 and anti--D5. 
The profile has the properties 



where m is an arbitrary integer. This picture, including the poles and periodicity, will 
be recovered with an improved quantitative description of r as a function of complex 
arguments in section 4 (see Figure 2). 

3.3 Solutions for the time-dependent fuzzy spheres: Collapsing 
spheres 

In the time-dependent case, there is a conserved energy 



r(a + 2E) 

r(a + 4m£) 
1 




r( (2m + 1)E ) 



0. 



(3.18) 



E 



(1 + r 4 ) 



(3.19) 



implying 




(3.20) 



This can be solved for the velocity 




(2r 4 + rg) - (2r 4 + r 8 ) 
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Writing s = ^ 



(1 + r 4 ) 2 - (1 + r 4 ) 2 

s = (i + r y < 3 - 2i; 

-{r A - rp{r A - rf) 
(1+r 4 ) 2 

A trivial redefinition s — > is relates the equation ()3.22j) to (J3.14)) . The time evolved can 
be written in terms of the radial distance 

f dt = I" dr - (1 + r ° 4) . (3.23) 

Jo Jro ^(2r 4 + rg)-(2r 4 + r8) 

Differentiating (j3.21j) gives a second order equation 

-4r 3 (l + r 4 ) 



E 2 

-4r 3 (l + r 4 



, 1 ■ (3-24) 

(1 + r 4 ) 2 

We can see from (j3.23|) that the time taken to start from r = r and reach r = is finite. 
We will call this finite time interval the time of collapse T. As in the spatial problem 
there are symmetries of the equation 

It : r(t) -> r(-t) 
I r : r(t) -> -r(t) 
T T : r(t) -> r(t - T) . (3.25) 

Following the same steps as in the spatial case, we can act successively with I^It^r , 
hulr, huh and then with integer multiples of 4T to produce a periodic solution defined 
for positive and negative time. As discussed before and in analogy to the Dl-D?> scenario, 
this can be interpreted in terms of collapsing-expanding ZM-branes. The radius as a 
function of t has the properties 

r(t + 4mT) = r(t) 
r(T) = r( (2m + 1)T ) = 

r{t + 2mT) = -r{t) , (3.26) 

where m is an arbitrary integer. We will see in the following that it will be useful to define 
a complex variable u\ whose real part is related to a and whose imaginary part is related 
to t as U\ = a — it. Unlike the case of the fuzzy two-sphere, we are not dealing simply 
with a variable u living on a torus. Rather it will become necessary to introduce a second 
complex variable u<i such that the pair (ui, U2) lives on the Jacobian of a genus two-curve. 
It will also be natural to impose a constraint which amounts to looking at sub- varieties of 
the Jacobian. At the end of this it will, nevertheless, be possible to recover the sequence 
of zeroes and poles in ()3.26|) and (|3.18j) . 
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3.4 Geometry and Automorphisms of Hyper-elliptic curve for 
the fuzzy S 

The equation ()3.22j) defines a Riemann surface of genus 3. The integrals of interest f)3. 
and ()3.23|) . are integrals of a holomorphic differential along certain cycles of the Riemann 
surface. It is useful to recall the Riemann-Hurwitz formula 



(2g - 2) = n{2G - 2) + B , (3.27) 

which gives the genus g of the covering surface in terms of the genus G of the target and 
the number of branch points B. Since the RHS of ()3.22|) is a polynomial of degree 8 there 
are 8 points where s = 0, i.e. 8 branch points. Here G = since the r co-ordinate can be 
viewed as living on the sphere, n = 2 and B = 8. So the integrals are defined on a genus 
g = 3 curve, which we will call S3. 

After dividing out by Tq on both sides of ()3.22|) and rescaling s = jf^i, r A = r 4 f 4 , 
rf = r 4 ff we have 

~ s 2 = (f 4 - l)(f 4 - ff) . (3.28) 
There are three independent holomorphic differentials on this curve coi = ^ ^2 = -f-, ^3 = 
(see for example P^ll5j). The infinitesimal time elapsed is dt = * = y . It will 
be important in calculating the integrals (j3.23j) ()3.16|) . to understand the automorphisms 
of the Riemann surface S3. In the limit of large tq, ff approaches —1 and there is an 
automorphism 




which leaves the equation of the curve unchanged. It transforms the holomorphic differ- 
entials as follows 



U\ — > iuj 3 
CU3 — > 



For any finite r there is a Z 2 automorphism, f — > —f. Quotienting by this can be 
achieved by changing variables f 2 = x. The first and third holomorphic differentials uj\ 
and CU3 transform as follows 

df dx 



^J(f 4 - l)(f 4 -ff) ^Ax(x 2 - l){x 2 -ff) 

f2df = Xdx . (3.29) 

V(f 4 - l)(f 4 -ff) v/4x(x 2 - l)(x 2 -ff) 
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We can view this in terms of a map from a genus three curve (r, s) to a genus two curve 



x = r 2 



y = 2fS, (3.30) 

which implies 

y 2 = 4f 2 (f 4 -l)(f 4 -f 4 ) = 4x(x 2 -l)(x 2 -f 4 ) . (3.31) 

This genus 2 curve (x,y) will be called E2. The holomorphic differentials are related by 
Y = — and ^j 11 = The map (j3.30|) has branching number B = 0, in agreement with 
()3.27|) . Even though || = at f = 0, this is not a branch point, since x is not a good 
local co-ordinate for the curve (x, y) at x = 0. Rather, a good local co-ordinate is y which 
is linearly related to f and this means that there is no branch point. 
The second holomorphic differential transforms as 



rdr dx 



^/(f 4 -l)(f 4 -rf) 2v/(x 2 - l)(x 2 -f 4 ) 
This can be viewed in terms of a map 



(3.32) 



x = r 2 



to a target torus (x,y), obeying 



y = s (3.33) 



y* = (x 2 - l)(x 2 - ?f) . (3.34) 



This map has B = 4, again in agreement with ()3.27j) . There are two branch points 
at (r,s) = (0, ±f 2 ) corresponding to (x,y) = (0, ±f 2 ). Similarly we have two branch 
points at f = 00. The region near f = 00 is best studied by defining variables f = f _1 
and s = sr~ 4 which re-express (I3.28J) as s 2 = (1 — f 4 )(l — f 4 f 4 ). Similarly we define 
x = = yx~ 2 which re-expresses ()3.34j) as y 2 = (1 — x 2 )(l — x 2 rf). Two branch 

points corresponding to f = 00 are at (x,y) = (0, ±1). 

The genus two curve £2 itself can be related to genus one curves. We first rewrite the 
(x, y) equation ()3.31|) as 

y 2 = Ax(x - l)(x + l)(x - iR\){x + iR\) , (3.35) 

with R\ = —if 2 and real. This is of the special form [TSl IT7] 

y 2 = x(x — l)(x — a)(x — j3)(x — a/3) , (3.36) 
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with a = — 1 and (3 = iR\. Such curves have an automorphism Ti 

X X 

3 



T l( ,) S y=<^ = ^, (3.37) 

X 6 X 6 



which transforms the holomorphic differentials 

,' dx\ iXdX 
Tx — 



/ xdx\ -i^fJidX 

Tl {—) = ~^r~ ■ (3 - 38) 

The genus two curve (|3.35|) can be put in the form ()3.36|) in yet another way, by choosing 
a = — 1, (3 — —iRf. This gives another automorphism T 2 which acts as 

?7? 2 



x 

3*7T 3 



T M = i^»L . (3.39) 

With either choice of (a, /3) one is led to look for variables which are invariant under the 
automorphism and as a result one describes the genus 2 curve as a covering of two genus 
one curves (£±,i]±) [To] IT7] 

4 = e ± (i-e ± )(i-4e ± ) with 4 = - ( ^Zf ] * - (3.40) 

(1 - a)(l - 0; 



The maps are given by 



(l-a)(l-/3)z 
(a; — a) (a; — /9) 

-V(l-«)(l-fl (j _'^gfl, »- (3.4D 



This gives an isomorphism of the Jacobian of the genus 2 curve in terms of a product of 
the Jacobians of the genus 1 curves £±. For the case a = — 1 we have = K'(k_). 

Since the complex structure is r = i^j^j, this means that the complex structures of S + 

and £_ are related by the SX(2, Z) transformation r — ► — -. Hence Ej- have isomorphic 
complex structures. 

As an aside we describe the full group of automorphisms of S 2 . It includes a, the 
hyperelliptic involution, a(x) = x, a(y) = —y ; T 3 , which acts as T^(x) = —x, T 3 (y) = —iy, 
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and T 4 , which acts as T 4 (x) = — x, T 4 (y) = iy. Relations in this group of automorphisms 
are 

TiT 2 = T 3 
T 2 T X = T 4 
T 3 T 4 = T 4 T 3 = 1 
T| = a 

T 3 a = T 4 . (3.42) 

In the limit r — > oo, i?i = 1 and the equation for the curve simplifies 

y 2 = x ( x 4 - l) . (3.43) 

As a result, the automorphism group is larger than at finite r . The automorphism group 
is generated by U\ and U 2 which act as follows 

1 



U x (x) 



Ui(y) 



x 

e 2 y 
x 3 



U 2 (x) = e 2 x 

U 2 (y) = e^y. (3.44) 

If we write U\(x) = X,Ui(y) = Y, we have the following action on the holomorphic 
differentials 

dx . XdX 

y y 

xdx ,dX 

= • 3.45 

y Y 

The action of U 2 on the holomorphic differentials is just 

dx e~<ix 



y y 

xdx e~^xdx 



(3.46) 

y y 

The automorphism group includes, as usual, the hyperelliptic involution acting as a(x) = 
x,a(y) = —y. There is also an element U 3 acting as U 3 (x) = —x,U 3 (y) = iy. There are 
relations Uf = a, U 2 = U 3 . In the large Tq limit, R x — > 1, the formulae for T X ,T 2 from 
(|3.37|) . ()3.39|) simplify and they can be written in terms of Ui,U 2 . Indeed we find that 

T x = aU 2 U x 
T 2 = aU 2 x U x . 
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3.5 Evaluation of integrals 

The time integral ()3.23|) can be done in terms of Appell functions. The indefinite integral 
is 

(l + r 4 )r f 1 1 1 5 r 4 r 4 



V(2r 4 + rg) Fl U ; 2'2 ; 4 ; r 4 ' 2 + r 4 J- (3 ' 47) 

A quick way to get this is by using the Integrator [THj but we outline a derivation. Ex- 
panding the integrand of (jHUHj) 



cir(r 4 -r 4 )-2(r 4 + 2 + r 4 )-5 



4fc / ^ 4 



^ r(f - *)r(* + i) r(f - z)r(/ + 1) v ; V2 + r, 

- a , 9 . 4,-1 ^ i r(i) rt|) 1 1 ,^r^ 4 V r 4 



fcmr(|-A;)r(i-Z)4A; + 4Z + l v ' W V2 + r 4 
Now use the following facts about T functions 

■1) 



ni-k) y ' r(|) 



■iy 



r(|-0 v ; r(j) 
A(k + i + \) r(f)r(§ + A + 

to recognise the series expansion of the Appell function [T9"] 



oo oo 



Fl (a; bl , b 2] c; z 1; = E E (a)fc+ M l)fc lm /Zl ' 2 > ( 3 - 48 ) 

with arguments as given in (J3.47)) and where we have also made use of the Pochhammer 
symbols (a) n = r{ ^ a f ■ 

Using ()3.47j) . the time taken to collapse from the initial radius r to a smaller radius 

r is 

(l + r 4 )r fl 1 1 5 r 4 \ (1 + r 4 )r fl 1 1 5 r 4 r 4 

v/(2r 4 + r«) 1 U'2'2'4' ' 2 + r 4 ^ > /(2r 4 + rg) 1 U' 2' 2' 4' r 4 ' 2 + r 4 , 

(3.49) 

For the special values r = r , Zi = 1, Ft of (|3.47j) simplifies to 



(1 + r 4 ) A 1 5 \ /l 13 



r v / (2T^j 2 i V4'2'4'V ^V4'2'4' 2 + r, 
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For r = 0, F\ — 1 and the indefinite integral ()3.47|) evaluates to zero. Hence the time 
of collapse is given by (|3.50|) which can be simplified to 

„_ (l + rg r(|)r(i) 13 -rj \ 

For large ro the last argument of the hypergeometric function simplifies to —1 and we get 

rp r(j)r(|) 2 1 4 r(|) r(|) 1iqQi3 / Q e<v\ 

It is also of interest to compute the interval in a along the D-string from the minimum 
size of the funnel cross-section to the place where the funnel blows up. The indefinite 
integral (|3.16|) gives by direct evaluation, just as above 

*(l + r 4 ) r (I 1 1 5 r 4 r 4 \ 

r 2 /tt—i 1 I 7' o' o' 7' H> "TTTTi ) • 



2 + r* V4 2 2 4 2 + 



The distance to blow-up is given by the difference of the last expression evaluated at 
infinity and at r . We do not have an exact formula for the large r asymptotics of the 
Appell function at finite r . We will thus be forced to take the large ro limit immediately. 
This will reduce (I3.53j) to 

(\ 1 9 r 8 \ 

and the distance to blow-up will be 

• r r(!)r(!) ir y/*rm 



A(-i) 1/8 r(f) 



(3.55) 



The final result is 



r Q (V2 - 1) V^fy ~ 0.4819..r . (3.56) 



Hence the full period 4E is 1.9276. .ro- The time of collapse is 1.1636. ..ro- The space and 
time periods are no longer the same as was the case for fuzzy S 2 , since there is a relative 
factor of (y/2- 1). 

4 Reduction of the g = 3 curve and inversion of the 
hyper-elliptic integral for the fuzzy-S* 4 

Whereas we have formulae for the time elapsed t in terms of r (13.49(1 or the Dl-co-ordinate 
a as a function of r by using ()3.53() . it is desirable to have the inverse formulae expressing 
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r as a function of a and t. It will turn out that, as in the case of the fuzzy S 2 , it will 
be useful to define a complex variable U\ = a — it. Whereas the u variable in the case of 
fuzzy S 2 lives on a genus one curve, here the story will involve higher genus curves and 
will require the introduction of a second complex variable w 2 . 

Let us revisit the cr-integral in the case of the fuzzy- S* 4 . The integral we want to 
perform is 

(1 + r 4 )<ir 



a 



with r\ = — (2 + r 4 ). If we make the re-scaling r = ^ and fi = we get 

0Tq r dr 



(l + ^o) Ji V(f 4 -l)(f 4 -f 4 ) 



(4.1) 



The RHS is the integral over the holomorphic differential on S 3 of section 3.4. Using the 
reduction to £ 2 by making a change of variables r 2 = x, we arrive at 



2tq a f x dx 



(1 + ^) Ji ^x[x 2 -\)[x 2 -ff) 
Similar steps for the time-dependent fuzzy-sphere give 



(4.2) 



= r dx (4 s) 

(1 + 4) Ji ^x(x 2 -l)(x 2 -rf) ' 

At this point it is useful to introduce a complex variable u\ = o — it. The inversion of 
the integrals ()4.2|) . ()4.3|) is related to the Jacobi inversion problem 

Xl dx f X2 dx 

— + — = Ml 

xo y j xo y 

Xl xdx f X2 xdx , , „. 
+ / = u 2 , (4.4) 

xo y j xo y 

where Xq is any fixed point on the Riemann surface E 2 . We will set Xq — 1. By further 
fixing X2 = 1 we recover the integral of interest in the first line 

Xl dx 

= Ml 

i y 

X1 T.d.T, 

— = m 2 , (4.5) 

'i y 

This is a constrained Jacobi inversion problem which is related to a sub-variety 4 of the 
Jacobian of S 2 , denoted as J(S 2 ). A naive attempt to consider the inversion of the first 

4 The geometry of such subvarieties is discussed extensively in j23 12] . One result is that H4.5JI defines 
a complex analytic homeomorphism from £2 to a complex analytic submanifold of J (£2)- 
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equation of f)4.5|) in isolation runs into difficulties with infinitesimal periods as explained 
on page 238 of [To] . 

By switching to the variables £±,r]± defined in (|3.41|) . the system (|4.4|) can be reduced 
to the sum of simple elliptic integrals 

+ / — = u+ 

u- , (4.6) 



6 [** 



where £i = £±(xi) and £ 2 = £±(^2)- We have used 

d£± /— r— -r, I — dx 



v/(l — a)(l — f3)(x ± a/o/3) — and u± = v/(l — a)(l — f3)(u 2 ± \fcif} u\) 

v± y 



The first of the two integrals can be brought into the form 

*^ 2dz 2dz 



h y/(l-Z 2 ){l-klz*) Jl x/{l-Z*){l-k 2 + Z 2 ) 

by the substitution £ = z 2 and then split to give 

^ dz dz 

+ 2 



(4.7) 
u + , (4.8) 



o v/(l-z 2 )(l-fc2* 2 ) 7 ^(l-^)(i_ fc 2 z2) 

_ 4 / ^ _ = u + , (4.9) 

7o v /(l-z2 )(1 _ A .2^ 2) 

which are just 

2Sn-\ x /l 1 , k + ) + 2Sn~\x/l 2 , k+) - AK(k + ) = u + . (4.10) 
Then by using the addition formulae for Sn^ 1 functions [22], we arrive at 



1 — k+£i« / 2 

(4.11) 

Thus by setting i 2 = 1 we get £2 = 1 and 

1 - 6 n 2 / M 



<?n 2 (^± + 2i^ + ),fc + ) . (4.12) 



1 - fc* & V 2 

After using the fact that Sn is anti-periodic in 2K(k) and implementing half-argument 
formulae, we end up with 

1 + Cn(u + ,k + ) = ^ 
1 + Dn(u + , k + ) 
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Starting from the second integral of ()4.6|) the same steps lead to 

1 + Cn(«_,L) 



1 + Dn(u-, k- 



6 • (4-14) 



These expressions hold the answer to the Jacobi inversion problem, although we still need 
to decouple u\ and u 2 in the u±'s. 

4.1 Series expansion of as a function of u\ 

The equations (j4.13j) and (j4.14j) imply that u + and w_ are constrained by 

= l + Cn(u + ,k + ) = 1 + Cn(u_,fc_) 
" 1 + Dn(u+, k+) 1 + Dn(u_,k_)' { ' ) 

This can be used to solve for 1*2(^1) (or tti^))- We do not have an explicit general 
solution to this transcendental constraint but we can solve it in a series around the initial 
radius x — 1. This corresponds to small times, i.e. corresponds to doing perturbation 
theory for u\ around zero and similarly for u 2 , as can be seen from the integrals (|4.5jl . 
We find 

U2 = u 1 + ^(l + Ri)ul + ^(l + Rf)(7 + SRt) U l 

+ 2^20 (1 + i? ' )(65 + + 9jR ' )M ' + 0{u " l] ■ (4 ' 16) 

As a consistency check, we can invert the two integrals independently and see whether 
the expansion ()4.16|) agrees with the results. 
In more detail we have 

dx 

1 — = = Ui , 

y/ X {x 2 -l){x 2 + R\) 

which can be expanded and inverted to give 

xi = \{l + R\)u\ + ^(1 + Rf)(7 + m\)u\ + ^L(l + J^)(65 + 66itf + 9i??K 

+-^—(1 + i^)(4645 + 7479/^ + 3087/2? + 189i# 2 )u? + 0(uf) . (4.17) 
The same can be done for 

xdx 

u 2 , 



1 y/x(x 2 -l)(x 2 + Rj) 
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giving 



x x = l(l + Ri)ul-^(l + R^(Rt-3)ui + ^(l + I^)(9 + 6Rt + 5Rl)ut 

+-^—(1 + i^)(189 + 63RI - 297RI - 235i?J 2 )^ + <D(uf) . (4.18) 

The expressions f!4.17|) and ()4.18|) can be combined to give the desired expansion of u 2 (u 1 ). 
We find perfect agreement between this and (j4.16|) . 

4.2 Evaluation of the time of collapse and the distance to blow- 
up 

We can use the previous discussion to give a new calculation of the time of collapse and 
the distance to blow-up, which can be shown to satisfy some non-trivial checks against 
the formulae in section 3. Consider (j4.15|) at the limits X\ = 1 and X\ = 0, where £i = 1 
and £i = 0. These lead to the two equations 

Cn(u±,k±) = and Cn(u±,k±) = -1 , (4.19) 

which give 

u x ± 1=1 = 2(2m± + l)K(k±) + Ain±K'(k±) and u x ± 1=0 = . (4.20) 

We should recall that for our case of a = —1, it turns out that k + is the complementary 
modulus of i.e. h^ + k^L = 1, which in turn means that 

K'(k ± ) = K(k T ) , (4.21) 

which is very useful in simplifying ()4.2()j) . From ()4.7|) we can write 

Ul = — + ^— - . 4.22 
2 v /(l-a)(l-/5) v / ^ 

Then by collecting like terms, the time of collapse or the distance to blow-up, will be 
extracted from imaginary or real values of U\ 

U\ = < 1=1 -< =0 

(2m + - 2m_ + l)K(k+) - (2m_ - 2in+ + l)K{kJ) 

y/-2(l-/3)P 

In order to compare this with something that we already know, we will consider the 
large-ro limit, where (3 — % and k\ = 1=F 2 ■ For these values we find that 

K(k-) = (V2 - i)K(k+) . (4.24) 
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Note that we should be careful here and observe a subtlety: By taking ro — > oo, we 
make k_ exactly real and larger than 1, while at very large values of ro, k_ has a small but 
non-zero positive imaginary part. It will be useful to recall the properties of the complete 
elliptic integral of the first kind K(k^), with k_ real. This function has a branch cut 
extending from 1 to oo. This means that the value above and below the cut will differ by 
a jump. This statement translates to [TH] 

lim K(k_ - ie) = K(k_) 

lim K(k_ + ie) = K(k) + 2iK'(k) . (4.25) 

Note that we are using the conventions for K(k±) of [22j[23 which differ from those of 
[T9] by a k — > k 2 . We can numerically verify that what we have corresponds to the second 
case and we will thus amend the above equation (|4.24jl to 

K{kJ) = (2i + V2-i)K(k + ) . (4.26) 

We are now ready to proceed. After some minor algebra and with the use of the tan ^ = 
1 + v^2, we reach 

Ul = -2- 3 / 4 cos - J (A + iB)K(k+) , (4.27) 
8 

where A = A\ + V2A 2 and B = B x + V2B 2 with 

At = 2m + + 2n_ + 2m_ - 6n + + 2 , A 2 = 2n_ - 2n + 

Bi = 2m + - 2n_ - 6m_ - 2n+ - 2 , B 2 = 2m + - 2m_ . (4.28) 

Now recall from (|4.2jl that since t is real and for the large-ro limit t = —roiii/2i, we need 
Ui to be purely imaginary in order to reproduce the results that we already have. This 
can be simply done by setting A = which is satisfied by 

n + = n_ 

m+ + m_ = 2n + - 1 (4.29) 

and has a first simple solution if we choose m_ = — 1 and m + = n + = n_ = 0. 
Then f!4.27j) becomes 

S-7T 

Ul = -i 2- 3 / 4 cos — (4 + 2V2)K(k + ) (4.30) 
8 

and with the relation u\ — a — it the time of collapse for large-ro follows 

cos ^:2- 3 / 4 (4 + 2v / 2) 




2.3272..— = 1.1636..r , (4.31) 
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exactly what we got previously from the evaluation of the integral in terms of Appell 
functions. When the conditions ()4.29|) for the vanishing of A are satisfied, the expressions 
for Bi and B 2 simplify to 

B 1 = -4(2m_ + 1) B 2 = -2(2m_ + 2n + + 1) . (4.32) 

Setting n + = we have Bi and B 2 both proportional to (2m_ + 1). This means that 
the first collapse to zero is repeated after every interval of twice the initial collapse time. 
This behaviour was anticipated using continuity and the symmetries (|3.26j) . 

Similarly, we can set B = and recover a real expression which will correspond to the 
distance that it takes for the funnel which has a minimum cross-section of tq to grow to 
infinity. We first use the large r$ limit. The reality conditions are 

m + = m„ 

n++n_ = -(2m_ + l) (4.33) 
and a simple solution is the one with m_ = m + = n + = and n_ = —1. Then 




roK | \l cos^2- 3 /V2 



~ 0.9639.. y = 0.4819..r , (4.34) 

again in agreement with the previous results. Generally when the conditions for the 
vanishing of the imaginary part are satisfied, expressions for A\ and A 2 simplify 

A x = -8n + A 2 = -2(2m+ + 2n + + 1) 

Setting n + = gives a sequence of poles at distances proportional to 2m + + 1. This 
was anticipated from symmetry and continuity using (J3.17|) in section 3 and is shown 
in Figure 2. Note that the pattern of zeroes along the time- axis is the same as for the 
fuzzy S 2 . The pattern of poles along the space axis is also the same as for the fuzzy S 2 . 
The difference is that the time from maximum radius to zero in the time evolution is not 
identical to the distance from minimum radius to infinite radius as in the case of fuzzy 
S 2 . There is, nevertheless, a simple ratio of (y/2 — 1) at the large r limit. 

It is also important to stress that we can use the above solutions to study the time of 
collapse for finite r and then compare with what one gets from (|3.51|) . We have checked 
numerically, for many values of r between zero and oo, that the time given in (|3.51|) 
agrees to six decimal digits with the expression 



; i + r 4 ) (K(k+) + K(k.) ) 
2r 3 v/-2/3(l-/3) 
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Figure 2: Numerical plot of the Jacobi elliptic function solution for the static S" 4 -fuzzy 
funnel array and the collapsing S^-sphere at large r . 



which reduces to (14.30)1 at large r . Similarly for the distance to blow-up we have 

E= (1 + tj) (2i + l)K(k + )-K(k-) 
2r 3 y/-2/3(l-/3) 

We believe the same numerical agreement to hold for finite r . However, since the 
asymptotics of (|3.53|) for large r at finite r were not well available with our mathematical 
software, we cannot confirm this. 

So far we have considered the time of collapse where r = ( or the distance to blow-up 
in the spatial case ), but we can also consider t as a function of r for any finite r in terms 
of the Appell function ([3.47)1 . The inverse expression of r in terms of t, and more generally 
the complex variable U\ = o — it, is contained in ([4.12)1 and the constraint ([4.15)1 . 



4.3 Solution of the problem in terms of the 112 variable and large- 
small duality 

We have seen that inverting the integrals ([3.23)1 and ([3.16)1 requires the definition of a 
complex variable u\ = <j\—it\ whose real and imaginary parts are related to the space and 
time variables. In addition we have to introduce the second holomorphic differential, so 
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that there are 2 complex variables m, w 2 . These variables are defined in terms of integrals 
and thus are subject to identifications by the period lattice C of integrals around the a 
and 6-cycles of the genus 2 Riemann surface, £ 2 . They live on C 2 /C, the Jacobian of 
E 2 . Introduction of a second point a; 2 on the Riemann surface relates our problem to the 
standard Jacobi inversion problem. The constraint a; 2 — 1 restricts to a subvariety of the 
Jacobian. So far the constraint ()4.15|) has been viewed as determining w 2 in terms of u\. 
This allowed us to describe function of u\ in the neighbourhood of X\ = 1 and 

also to get new formulae for the time of collapse/distance to blow-up which have been 
checked numerically against evaluation of the integrals using Appell functions. However, 
it is also of interest to consider using the constraint to solve U\ in terms of w 2 and hence 
describe function of w 2 . The reason for this is that the automorphisms of the 

Riemann surface allow us to relate the large r (equivalently large X\) behaviour of the 
spatial problem described in terms of the u\ variable, to the small r (equivalently small 
X\) behaviour of the time dependent problem described in terms of the w 2 variable and 
vice versa. The relation is simpler in the large r limit, hence we describe this first and 
then we return to the case of finite r . 

The use of the w 2 variable is natural if we introduce another Lagrangian of the type 
(|B.1[) . with a = 2. The original Lagrangian of interest with a = and the a = 2 
Lagrangian are coupled in the Jacobi Inversion problem f)4.4j) as well as the constrained 
version of the Jacobi Inversion problem (J4.5)) obtained by setting x 2 = 1. Just as u\ is the 
complexified variable for the first Lagrangian, w 2 is the complexified space-time variable 
for the second Lagrangian. It will be convenient, for the following discussion, to define 
the relation between the second set of space-time variables t 2 , <r 2 and the complex variable 
u 2 as u 2 = — cr 2 + it 2 . 

We have also seen that the automorphism U\ of section 3.4 maps du\ to —idu2 as 
in ()3.45j) . hence takes the time dependent problem for the first kind of Lagrangian to 
the static one for the second mapping zeroes of a given periodicity to poles of the same 
periodicity. We will check this by investigating what happens to w 2 when we follow the 
collapse of x\ down to zero along imaginary w 2 , or the blow-up of x\ to infinity along real 
m 2 . It is clear from the integrals in (|4.5|) (and the series expansion (|4.16Jl ) that when x\ 
decreases from 1 to 0, both U\ and w 2 are imaginary and when x\ increases from 1 to 
infinity, both are real. 

Solving for m 2 , the second of the equations ()4.7|) will give 



u + + U- 



(4.37) 



u 2 = 



2y/(l -a)(l-P) 



and by following the steps leading to 1)4.27)1 and (J4.28)) we arrive at 



u 2 = - 



2- 3 / 4 cos --{A 1 + iB')K(k+) . 



(4.38) 
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Here we used £an| = \/2 — 1. The A's and B's are given by 

A\ = 2m + + 2n_ + 2m_ - 6n + + 2 , A' 2 = 2n + - 2n_ 

B[ = 2 + 6m„ + 2n_ - 2m + + 2n + , 5 2 = 2m + - 2m_ . (4.39) 

Notice by comparing P^Kjl and (jOHjl . that = A 1 ,B' 1 = -B x and A' 2 = -A 2 , B' 2 = B 2 . 
By choosing m_ = — 1, m + = n_ = n + = we get 

m 2 = 2z(2 - V2) cos ^2- 3 / 4 ^ (fc + ) 

8 

~ i 0.9639.. , (4.40) 

which is exactly the same as in (|4.34|) in agreement with what we expect from the auto- 
morphism. Similarly for the real case 

u 2 = -2v / 2cos^2~ 3/4 ir(A; + ) 
8 

~ -2.3271.. . (4.41) 

With the definition u 2 = —cr 2 + it 2 we see that we are getting the expected positive time 
of collapse in terms of the alternative time-variable and the expected positive distance 
to blow-up with the alternative spatial variable. The matching of the A' and B 1 with 
A and B guarantees that x±, expressed as a function of u 2 , will have zeroes along the 
imaginary axis at odd integer multiples of a basic time of collapse and poles along the 
real axis at odd integer multiples of a basic distance to blow-up. The time of collapse for 
u 2 is the same as the distance to blow-up for u\ and the distance to blow-up for u 2 is 
the same as the time of collapse for u\. This gives a precise map between the behaviour 
at zero and the behaviour at infinity of the radius of the fuzzy sphere, generalising the 
relations that were found for the case of the fuzzy 2-sphere. This relation is expected 
from the automorphism JJ\ for large Tq described in section 3, which maps U\ to —iu 2 
f)3.45p . This precise large-small relation for space and time-dependent fuzzy spheres is 
physically interesting. The physics of the large r limit is very well understood because 
it corresponds to the D-strings blowing up into a D5-brane. The physics of the small r 
limit appears mysterious because it involves sub-stringy distances. We have shown that 
the two regions are closely related through the underlying Riemann surface which unifies 
the space and time aspects of the problem. 

4.3.1 Large-small duality at finite r 

We saw in the discussion above that the spatial problem of the fuzzy S 4 evolving from 
a minimum size at r to infinity can be related to the time-dependent problem of the 
fuzzy sphere collapsing from r to zero. This large-small relation involves a map between 
the u\ description and the u 2 description of the problem and uses a simplification which 
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is valid at large r . There continues to be a large-small duality at finite r , but it is 
slightly more involved than the one at large r . The difference is due to the nature of the 
automorphisms of the Riemann surface at finite r and in the large Tq limit, which were 
described in detail in section 3.4. Indeed, the discussion above used the automorphism 
U\ in a crucial way. 

To describe the duality at finite Tq, it is useful to think of a problem similar to the 
one we considered above, by choosing a different basepoint in ()4.5|) . i.e. x = —iR\ 

El dx 
ml V 

= u 2 . (4.42) 



iR\ y 

The upper limit is chosen, in the first instance, to vary along the negative imaginary axis 
up to zero. A second problem is to consider x\ extending along the negative imaginary 
axis down to infinity. 

Applying the automorphism T\ of ()3.37|) to the first line of ()4.5|) we have 



Mi(Xi) 



X1 dx 

1 y 

e f /■*!=— s 1 XdX 



Ri J-iRi y 

/ iR\ 
-u 2 xi - 



Similarly we have 



Ri \ x 



iR 2 ' 

u 2 (x 1 ) = Rie^ui ( xi = ) . (4.43) 



As we saw earlier in this section, the solution to (|4.5jl can be described in terms of either 
the u\ or the u 2 variable. Likewise the inversion of 1)4.42)1 can be expressed in terms of 
either u\ or u 2 . The action of the automorphism T\ described above implies that the 
solution of the spatial problem (J4.5)) where x\ evolves from 1 to infinity along the real 
axis, when given in terms of the real part of the u\ variable, maps to the evolution of x\ 
from — %R\ along the imaginary axis to zero, as described by the u 2 variable. Similarly 
the time-dependent problem of x\ evolving from 1 to zero when described in terms of the 
imaginary part of the u\ variable, maps to the evolution of X\ along the imaginary axis 
from —iRl to infinity, as described by the u 2 variable. This shows that there continues 
to be a large-small duality at finite ro, but it relates the original problem with real x to 
a problem with imaginary x. 
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5 Space and Time dependent Fuzzy-/S 6 

Here we will briefly talk about the nature of the solution in the case of the Dl _L D7 
intersection, discussed by [3], which involves the fuzzy-S* 6 . Starting from the D-string 
theory point of view, we will now have seven transverse scalar fields, given by the ansatz 

$*(<7,t) = ±R(a,T)G\ z = 1,...7, (5.1) 

where the G l, s are given by the action of the 5*0(7) P's on the symmetric and traceless 
n-fold tensor product of the basic spinor V, the dimension of which is related to n by 



(n + l)(n + 2)(n + 3) 2 (n + 4)(n + 5) 

360 ' 1 ' ' 

Again, the radial profile and the fuzzy-S* 6 physical radius are related by 

R P h(v, t) = y/c\R(a, r) , (5.3) 

with c the quadratic Casimir G % G l = cIn = n(n + 6)1nxn- The time-dependent general- 
isation for the leading 1/N action of [3] can be written in dimensionless variables, which 
are once more defined as in (|2.10jh 

S t = -NT x J d 2 aVl + r' 2 - f 2 (l + r 4 ) 3/2 . (5.4) 

In a manner identical to the discussion in section 3, the equations of motion can be given 
in a Lorentz-invariant expression 



1 + (<V) (#*r) 



dpdPr + (d^r) (d v r) (d u r) - (d^r) (d u r) (d^r) = 6 r A ^ \ + r l J ■ ( 5 - 5 ) 

At this point it is natural to propose that the form of the action and the equations of 
motion will generalise in a nice way for any fuzzy-5' 2fc sphere 



Si = -TxJ da 2 STr\/(l + r' 2 - r 2 ) (1 + r 4 ) K , (5.6) 
with the large- N equations of motion 



d^r + {d^r) (d u r) (d u r) - {d^r) (d v r) (0"r) = 2k r 3 ^ ^ ^ j • (5.7) 

As we saw in the previous cases, there will be a curve related to the blow-up of the 
funnel, derived by the conservation of pressure if we restrict to static configurations and 
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also to the corresponding collapse of a L>6-brane by conservation of energy if we completely 
drop the space variable. We find that the curve determining the solutions is 

s 2 = (r 4 - r 4 ) (r 4 - rf) (r 4 - rf) , (5.8) 

which is of genus 5 and where a factor of (1 + r 4 ) has been absorbed in the definition of 
s. The roots are given by 

(l + r 4 ) = 4 
(1+r 4 ) = 4v 

(1+r 4 ) = 4v\ (5.9) 



with rj = exp 

• Automorphism at large r : 

At large r , we have rf = r^rj and rf = r^rf. Then there exists an automorphism 



R 

S 



r 2 

r 

_> _ s2r l 2 



r 12 



It is convenient to define f — —, and s 2 = s 2 rl 2 . In these variables 

r ' 



1 

f 

is 



and the action on the holomorphic differentials is 

uji — > icj 5 
cj 2 — > ?CJ 4 

cj 4 — > icj 2 
u;5 — > ?cji . 



Automorphism at r = : 

Now we have 



= r 4( r 4_ r 4 )(r 4_ r 4 ) (51Q) 
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and there is an automorphism 



R riT2 
r 
.-2 



2 S 2 {r 1 r 2 f 



^.16 

In this limit rf = (77 — 1) and rf = (r? 2 — 1), so in the formulae above we can write 
r\r\ = 3. 

Symmetry at finite ro : 

It is useful to write the curve in terms of a variable u defined by u A = (1 + r 4 ) and 
to write w 4 . = 1 + r 4 . Then we have 

s 2 = (u 4 -ut)(u 4 -utv)(u*-uW). (5.11) 

2 2 12 

A symmetry is v = — , s 2 = — s "° . Expressing the symmetry in the r, s variables, 
we have i?, S 1 obeying the same equation (|5.8|) with 

(1 + * 4 ) = il + r » )2 



S 



(1 + r 4 ) 



1+r 



4A3 



This reduces to the R = ^ for r >> 1, i? >> 1, ro >> 1. Unfortunately R is not a 
rational function of r, but an algebraic function of r involving fourth roots, hence 
it is not a holomorphic or meromorphic function. Hence it is not possible to use 
this symmetry to map the holomorphic differentials of the genus 5 curve to those on 
genus 1 curves. We can still make the change of variables x = r 2 to get a reduction 
down to genus 3, but we have not been able to reduce this any further. 

Since we cannot reduce the curve down to a product of genus one curves, we cannot 
relate the problem of inverting the hyper-elliptic integral to elliptic functions. We can 
nevertheless relate it to the Jacobi inversion problem at genus 3. We consider variables 
ui,u 2 ,u 3 defined as 

UJi + / 0Ji+ UJx = Ul 

! Jl Jl 

XI PX2 PX3 

UJ 2 + UJ 2 + UJ 2 = U 2 

xi i-x 2 i-x 3 

LU 3 + J LU 3 + J LU 3 = U 3 . 
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The variables u\,u 2 ,u 3 live on the Jacobian of the genus three curve, which is a 
complex torus of the form C 3 /£. The integrands appearing above live naturally on a 
Riemann surface which is a cover of the complex plane, branched at 8 points. The lattice 
C arises from doing the integrals around the a and 6-cycles of the Riemann surface. The 
equations (|5.12|) can be inverted to express 



Xi = x 1 (u 1 ,u 2 , u 3 ) 
x 2 = x 2 (u 1} u 2 ,u 3 ) 
x 3 = x 3 (u 1 ,u 2 ,u 3 ) , 

where Xi(ui,u 2 ,Us) can be given in terms of genus three theta functions, or equivalently 
in terms of hyper-elliptic Kleinian functions j2S]- The system (|5.12|) simplifies if we set 
x 2 — x 3 — 1. These simplified equations define a sub- variety of the Jacobian which 
is isomorphic to the genus 3 Riemann surface we started with [20] ■ The constraints 
x 2 = x% = 1 can be used to solve, at least locally near x\ = 1, for u 2 , u 3 in terms of u\. 
Then we can write xi(ui,u 2 ,u 3 ) as X\{ui). This program was carried out explicitly in 
section 3, where the higher genus theta functions degenerated into expressions in terms 
of ordinary elliptic functions thanks to the reduction of the genus three curve we started 
with, to a product of genus one curves. For completeness, we give a short review of the 
solution to the Jacobi Inversion problem in terms of higher genus $ functions and the 
related Kleinian cr-function. 



5.1 Jacobi inversion problem and d functions 

The general Jacobi inversion problem can be formulated as follows fHl EH I2S] ■ F° r an y 
hyper-elliptic curve E of genus g, realised as a 2-sheeted cover over a Riemann sphere 

2g+2 g+1 g+1 

i=0 i=l i=l 

with a's and c's being the branch points, between which we stretch the cuts, the system 
of integral equations 




describes the invertible Abel map, U : S s / 'S g — > Jac(E), taking g symmetric points from 
the Riemann surface to the Jacobian of E. The latter is just C 9 /C, where £ = 2( © 2£' is 
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the lattice that is generated by the non-degenerate periods of the holomorphic differentials, 
or differentials of the first kind, defined on the surface 



2(= f uji and 2(' = & 
JAv. Jb, 



with 



x % l dx 



UJi = 



y 



i = l,...,g . 



(5.13) 



(5.14) 



The period matrix is the g x g matrix given by r = C _1 C' an d belongs to the Siegel upper 
halfspace of degree g, having positive imaginary part and being symmetric. There is also 
a set of canonical meromorphic differentials, or differentials of the second kind, naturally 
defined on the Riemann surface 



2<7+l-j 
k=j 



x k dx 



3 = l,---,9, 



(5.15) 



the periods of which are 



2i] = - (b £, and 2rf = - <b & . 

J A k JBi. 



Now consider m, m' e Z s two arbitrary vectors and define the periods 

Q(m, to') = 2(m + 2('m' and E(m, m') = 2r]m + 2r\vn! . 



(5.16) 



(5.17) 



We can define the fundamental Kleinian a-function in terms of higher genus Riemann 
^-functions 

a(u) = Ce uTKU ${{2u- l )u- K a \r) , (5.18) 

where C is a constant, k = (2()~ 1 rj and K a is the vector of Riemann constants with base 
point a, given by the Riemann vanishing theorem 



k=i 



(5.19) 



with dv = (2( 1 )(c^i, • • • ,^g) T being the set of normalised canonical holomorphic differ- 
entials. Furthermore, the genus-g Riemann ^-function with half-integer characteristics 

[e] = 



' e' ' 






■ 4" 


e 




. £l ■ 


. Eg _ 



E C 29 is defined as 

7?[e](«|r) = ex P ni K m + £ '^ T ( m + £ ') + 2 ( v + £ ) T ( m + • ( 5 - 20 ) 
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The general pre-image of the Abel map can be given in a simple algebraic form as the 
roots of a polynomial in x, given by 



P(x; u) = x 9 - x 9 1 p g ,g{u) - x 9 2 p g , g _ 1 (u) - ... - p g>1 {u) . (5.21) 

Here u — (ui, . . . , u g ) T and the p's are higher genus versions of the standard Weierstrass 
elliptic p-functions, which are defined as the logarithmic derivatives of the fundamental 
hyper-elliptic cr-functions 

. . d 2 In a(u) . . <9 3 In a(u) . . , . nnS 

M«) = ~^-. ^ ) = "^^'---' «^*.-=l.-.ff ( 5 - 22 ) 

and have the nice periodicity properties 

fPij(u + Q(m, m')) = pij(u) , i,j = l,...,g. (5.23) 

We can make use of all this technology and give at least an implicit solution for the S* 6 at 
the stage where we have been able to reduce the problem from one of genus-5, to one of 
genus-3. By equating the general polynomial with roots Xi,x 2 , x 3 to the polynomial with 
^-coefficients, we get 

Xi + X 2 + X 3 = p 33 (Ul, «2, W3) 
Xl X 2 X 3 = p23(u 1 ,U 2 ,U 3 ) 

xix 2 + x 2 x 3 + = -pi 3 (ui, u 2 , u 3 ) . (5.24) 

As we have already mentioned, we can fix two of the three points to be x 2 — x 3 — 1, to 
get 

xi + 2 = p 33 {ui,u 2 ,u 3 ) 

Xl = p23{Ul,U 2 ,U 3 ) 

2xi + 2 = -pi3(u u u 2 ,u 3 ) . (5.25) 

This implies two transcendental constraints which can be used to get a solution for x±, in 
terms of the p's as functions of u±. 



6 Summary and Outlook 

The space and time dependence of fuzzy spheres S 2 , S* 4 and S 6 are governed by equations 
which follow from the DBI action of D-branes. These fuzzy spheres can arise as bound 
states of DO with D2, D4 or D6 on the one hand, or as cross-sections of fuzzy funnels 
formed by Dl expanding into D3, D5 or D7 on the other. The purely time-dependent pro- 
cess has the simplest realisation as the collapsing-expanding DO — _D2p-brane, although it 
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also arises as a time dependent process for the Dl — D(2p+ 1) equations with no variation 
along the spatial Z)l-direction. The first order equations of motion are closely related to 
some Riemann surfaces and the infinitesimal time elapsed or the infinitesimal distance 
along the Dl, are given by a holomorphic differential on these surfaces. We have given 
descriptions of the solutions in terms of elliptic functions or their higher genus gener- 
alisations. We showed that the space and time processes exhibit some very interesting 
large-small dualities closely related to the geometry of the Riemann surfaces. 

In the case of S 2 , the Riemann surface has genus one and the duality is related to 
properties of Jacobi elliptic functions known as "complex multiplication formulae" [Hj. The 
genus 1 Riemann surface that arises here has automorphisms, holomorphic maps to itself, 
which are responsible for these properties. We observed that the large-small duality is also 
directly related to a transformation property of the second order differential equations. 

In the case of S* 4 we derived formulae for the distance a(r) along the Dl-brane as a 
function of the radius of the funnel's fuzzy sphere cross-section and the time elapsed t(r) 
as a function of the radius. These formulae were expressed in terms of special functions 
such as Appell and hypergeometric functions. It was useful to combine space and time 
into a complex variable u\ = a — it which was related to a genus three Riemann surface, 
having a number of automorphisms. These automorphisms were used to eventually relate 
the problem to a product of two genus one Riemann surfaces. After introducing a second 
complex variable w 2 , the problem of inverting the integrals to get a formula for the radius 
as a function of the complexified space-time variable u\ was related to a classical problem 
in Riemann surfaces, the Jacobi Inversion problem. This can be solved, in general, using 
higher genus theta functions. The relation to a product of genus one surfaces in the case 
at hand, means that the inversion problem can be expressed in terms of the standard 
elliptic integrals. This allowed us to give a solution of the inversion problem in terms of 
standard Jacobi elliptic functions. The solution involves a constraint in terms of elliptic 
functions, which can be used to solve u<i in terms u±, or vice versa. This approach yields 
a new construction of r(ui) as a series around r = tq, which agrees with a direct series 
inversion of the Appell function. It also gives new formulae for the time of collapse and 
distance to blow-up in terms of sums of complete elliptic integrals, which were checked 
numerically to agree with the formulae in terms of hypergeometric functions. 

The automorphism which allows a reduction of the problem to one involving genus 
one surfaces also relates the large r behaviour of the spatial (time-dependent) solutions 
to the small r behaviour of the time-dependent (spatial) problem. The introduction of 
the extra variable U2, required to make a connection to the Jacobi inversion problem, also 
enters the description of this large-small duality. 

Our discussion in the S 6 case was less complete, but a lot of the structure uncovered 
above continues to apply. The integrals giving t(r) or a(r) are integrals of a holomorphic 
differential on a genus 5 Riemann surface. A simple R = r 2 transformation maps it to 
a holomorphic differential on a genus 3 curve. The inversion problem of expressing r in 
terms of u± = a — it can be related to the Jacobi inversion problem and a solution in terms 
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of higher genus theta functions is outlined. We did not find any automorphisms of the 
genus 3 Riemann surface which would relate the problem to one involving holomorphic 
differentials on a genus 1 curve. There do continue to exist symmetry transformations of 
the curve but they are no longer holomorphic in terms of the (r, s) variables. 

It will be interesting to consider the full CFT description of these funnels. Some steps 
towards the full CFT description of these systems has been made There should be 
a boundary state describing the spatial configuration of a Dl-branes forming a funnel 
which blows up into a D3-brane. One could start with the CFT of the multiple Dl- 
branes and consider a boundary perturbation describing the funnel which opens into a 
D3. Alternatively we could start with a CFT for the D3 and introduce boundary per- 
turbations corresponding to the magnetic field strength and the transverse scalar excited 
on the brane, which describe the Dl-spike. In either case the boundary perturbation 
will involve the elliptic functions which appeared in section 2. Similarly the boundary 
perturbation for the fuzzy S 4 case would involve elliptic functions associated with a pair 
of genus 1 Riemann surfaces, of the type described in section 4. One such boundary 
state corresponds to the purely time-dependent solution and another to the purely space- 
dependent solution. If the large-small duality of the DBI action continues to hold in the 
CFT context, this would give a remarkable relation between the zero radius limit of time- 
dependent brane collapse and the well-understood blow-up of a Dl-brane funnel into a 
higher dimensional brane. Analytic continuations from Minkowski to Euclidean space in 
the context of boundary states have proved useful in studies of the rolling tachyon |27j . 
We expect similar applications of the analytic structures described here. 

A complementary way to approach these dualities is to consider the supergravity de- 
scription. The time dependent system of a collapsing D0-D2-brane, for example, should 
correspond to a time-dependent supergravity background, when the backreaction of space- 
time to the stress tensor of the D0-D2 system is taken into account. A Wick rotation 
of this background can convert the time-like co-ordinate t to a space-like co-ordinate. 
Similarly the spatial funnel of a D1-D3 system has a supergravity description with an S 2 
having a size that depends on the distance along the direction of the D-string. The Wick 
rotated D2-background should be related by a large-small duality to the D1-D3 system. 
This is a sort of generalisation of S 1 T-duality to S 2 . Similarly, the DO — D4 system 
is related to the Dl _L D5. The Jacobi Inversion problem which underlies our solutions 
has already arisen in the context of general relativity [28J. It will be interesting to see if 
this type of GR application of the Jacobi Inversion also arises in the spacetime solutions 
corresponding to our brane configurations. Another perspective on the geometry of these 
configurations , including the boosted ones in the Appendix A.l, and on the relations 
to S-branes, should be provided by considering the induced metric on the brane as for 
example in [2*9"] . 

Generalisation of this discussion to funnels of CP 2 and other cross-sections, as well as 
to funnels made of dyonic strings will be interesting. The study of finite n corrections as 
in 0] is another interesting avenue. 
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Studying fluctuations around the solutions considered is likely to involve equations 
where the elliptic functions and their higher genus generalisations appear as potentials in 
Dirac or Klein-Gordon type equations for the fluctuations. Schrodinger equations with 
such potentials are well studied in the literature of integrable models 30J. It will be 
interesting to see if such integrable models have a role in the analysis of fluctuations. 

Finally, the genus 3 curve for the S 4 was very special. It is a related by the R = r 2 
map to a product of genus 1 and a genus 2 curve. The resulting genus 2 curve itself has a 
Jacobian which is related to K3 Kummer surfaces and factorises into a product of genus 
one curves [TH]- Special K3 J s and complex multiplication also appear in the attractor 
mechanism [31] E2]- Whether the appearances of special K3's in these two very different 
ways in string theory are related is an intriguing question. 
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A Further aspects of space-time dependence of S* 2 
funnel 

A.l Lorentz Invariance of the BPS condition 

We will study the supersymmetry of the space and time-dependent system considered in 
section 2, in order to obtain a Lorentz invariant BPS condition. The vanishing of the 
variation of the gaugino on the world-volume of the intersection will require that PP 

6 X = T^F, u e = 0, (A.l) 

where //, v spacetime indices, with 

F ab = F ai = D a ^ I], /<I>'.<K and a = a,r , z = l,2,3. (A.2) 

Here the spinor e already satisfies the D-string projection T Tcr e = e and multiply indexed 
T's denote their normalised, fully antisymmetric product, e.g. Y^ u = ^[T^,T U ]. Then 

(iT jk [&, + 2T ai D a &) e = . (A.3) 

We will consider the conjugate spinor 

^ = 5 X t r° (A.4) 
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and construct the invariant quantity 

5 X ] r° 5 X = , (A.5) 
where our conventions for the T-matrices are 

(r ) 2 = -l , (r) 2 = 1 and (r^) f = r°rT , (r^) f = r°rT° . (A.6) 

Then ()A.5|) becomes 

e t (pftR)t[$< j $™] + 2(r fei ) t J D f) $ l ) r° (i F fc [$ j , $ fe ] + 2r M J D a $ i ) e = . (A.7) 

This will give three kinds of terms: Quadratic terms in derivatives, a quadratic term in 
commutators of and cross-terms. The first of these will be 

4(D a ¥)(D b ¥) e t (r°r fe 'T )(r r ai )e . (A.8) 

By making use of the ansatz = R(a, r)a l and by taking a trace over the gauge indices 
this becomes 

A(d a h)(d b R) Tr (aV) e f (V i{r b , T a }T j r + r°^[r fe , r a ]F'r^ e (A.9) 

Since {T^, T u } = 2g tiu , Tr a l a^ = c 8 % i for a constant c and the expression is symmetric 
in a, b we get 

12 c(d a R)(d a R) (e+rV) . (A. 10) 

The product of commutators will be 

-[^,$ m p J ',$ fe ] e tr r £m r J \ , (A.ii) 

which with the trace and the use of the ansatz becomes 

4ci? 4 e tmp e jkq 5 pq ^T°T em T jk e . (A. 12) 

Since we are looking for a supersymmetry condition for the intersection, we want the 
spinor e to satisfy p" 1 - 7 e = e lJ e appropriate for the D3-brane. Using this fact we can 
re-write the above as 

-4ci? 4 {5 l] 5 mk - 5 ek 5 mj ) e imn e jkn e f r°e (A. 13) 

and end up with 

-12c(4R 4 ) (e^e) . (A. 14) 

Finally, the cross-terms 

2% (e j r°r jk r ai [<s> j , (£> a $> - e t r Q r ai r ifc (^ a $ i )[$ i , $ fc ] e ) (A.15) 
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will vanish by following the same steps, since Tr a l a^a k ~ e 1 ^ and the T-matrix commu- 
tator is [P,P' fc ] = (5 ij T k - 5 ik Ti). Collecting (IA~TU1) and (lA~T4j) we have 

e f r°e ((d a R)(d a R) - 4 i? 4 ) = . 

which implies 

(d a R)(d a R) -AR A = (A. 16) 

By employing dimensionless variables, the supersymmetry-preserving condition is simply 

(<V) (#V) = r 4 , fi = a,T. (A. 17) 

Thus, every configuration that satisfies ()A.17|) qualifies as a BPS state and therefore will 
be semi-classically stable. The special case R' = 2(R) 2 was used in JTj and is related to 
Nahm's equations. 

We will briefly look at some solutions of this first-order equation. It is not hard to see 
that the simplest ones can be put into the form 

2 7(cr - (Too + pt) 

where o"oo a constant, 7 = ^-3— _ and /3 another constant, which can be thought of as 

the ratio of the velocity of the system over the speed of light c—1. Then the solution at 
hand is nothing but a boost of the static funnel 

R = ± W - r- (A.19) 

2(cr - (Too) 

Moreover (|A.18J) will satisfy the full non-linear ()2.9)1 and is therefore a proper solution of 
our space and time-dependent system. In this case we can see explicitly how this solution 
should also satisfy the YM equations of motion. By starting with the BPS-like relation 
(IA.17|) . one gets by differentiation and algebraic manipulation 

dfJSTr (d v r)(d^r) = 2r 3 (d II r)(d> M r) , (A.20) 

which if substituted in (|2.12j) and again with (jA.17|) . will yield exactly d^d^r = 2r 3 , which 
are the YM equations. Relations between the BPS condition, the DBI and YM equations 
were discussed in jHl]. Therefore the class of space and time- dependent supersymmetric 
solutions to DBI simultaneously satisfy the BPS and the YM equations. Interestingly 
there exist solutions to one of YM or BPS which do not solve the DBI. For example, 
r Y M = / 1 solves the YM equations but not the BPS or DBI and t BP q = -i=k — * is 

just a solution to BPS. 

The fact that we have recovered, in the context of space and time dependent trans- 
verse scalars, a boosted static solution is something that we should have expected: The 
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expression for the Born-Infeld action is, of course, Lorentz invariant. This guarantees that 
its extrema, and solutions to the equations of motion, although transforming accordingly 
under boosts, should still be valid solutions at the new points x' = A~ 1 x. The lowest 
energy configuration of the system is naturally the BPS one and a boost provides a gen- 
eralisation that is still stable and time-dependent. A natural consequence of this is the 
boosted brane array, which is indeed a solution of the space and time dependent DBI 
equations of motion 

r(a, r) = ±r ^ (A.21) 



Cn 

and the boosted collapsing brane 



V2r -y(a- (3t) 1 




A. 2 Chern-Simons terms and D-brane charges for D1-D3 sys- 
tem 

We found in section 2.1 that the equations of motion for a space and time-dependent 
'funnel' configuration have a solution both in the Dl and the D3 picture. The Chern- 
Simons part of the non-abelian DBI in a flat background, is not relevant in the search 
for the equations of motion. However, it would be a nice check to see whether the charge 
calculation also agrees when time dependence of the scalars is introduced, as happens for 
the static case. We will be assuming spherical symmetry of the solutions in what follows. 

We expect to recover a coupling to a higher dimensional brane through the non- 
commutative, transverse scalars and the dielectric effect. This is indeed the case and 
one gets for the D3-charge of the Dl configuration for any time-dependent solution, at 
large- 

= ^JdrAda R\ h (C^k ph - C { XKh) ■ (A.23) 
The Chern-Simons part of the low energy effective D3-brane action reads 

S C s = ^J Str P[C {4) ] + fi 3 J Str P[C {2) ] A F + fi 3 J Str P[C {0) ] A F A F . (A.24) 

By considering a spherical co-ordinate embedding in the static gauge, the calculation of 
the D3 charge will give 

Sgl = 4vr/i3 JdtA dR D3 R 2 m [c^ + Cg^A*] , (A.25) 

where Rd3, t are world- volume indices and 1, 2, 3 are space-time indices. 
Following [2j we will make the following identifications 
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• The physical radius of the fuzzy S 2 on the Dl side should correspond to the co- 
ordinate i?£>3 on the world- volume of the D3 

R P h <— > Rd3 • (A. 26) 

• The a co-ordinate on the world- volume of the Dl should be analogous to the trans- 
verse scalar $ on the D3 

o < — ► A$ . (A.27) 

• We finally assume 

r < — ► t . (A.28) 

By thinking about the physical radius of the fuzzy sphere as a function of a and t, 
using the relationships (|2.18j) and also the fact {2 r n\fa'Y~ p fjP' = fjP, we are able to write 
(jA.25|) in terms of Dl world- volume quantities. Everything works out nicely and one gets 

SSk = SE! ■ (A.29) 
Furthermore we can make use of (|A~T8|) to get: 

$(7W) = ± (V^ - ^rl ( A - 3 °) 



as the boosted solution on the D3. This can be re-written as 

A$(i? D3 , t) = 7 (A$ (i?D3) + 0) , (A.31) 

where t' is related to t by a Lorentz transformation t' — "fit — /3A$) and can be viewed as 
the time co-ordinate for an observer on the worldvolume of the boosted D3. 



B Lagrangians for holomorphic differentials 

We have seen that the space or time-dependence of the radius of the fuzzy S 4 described 
by the Lagrangian (|3.4|) is given by integrating the holomorphic differential — on the 
curve (r, s) given by ()3.22|) . There are more general holomorphic differentials on the same 
Riemann surface, which enter on equal footing in the geometry of the Riemann surface. 
For the hyperelliptic curves of the type we considered, they are of the form = dt a , 
where a can take values from 1 to the genus of the curve. It is natural to ask if there 
are Lagrangians such that the time elapsed or distance are given by the more general 
holomorphic differentials. This is indeed possible and the Lagrangian densities are 

L a = -sjl -r 2 «(r 2 + r /2 )(l + r 4 ) . (B.l) 
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The automorphisms we discussed in section 4, can be used to relate the time evolution 
of interest which appears as the imaginary part of ui, to the space evolution for a = 2. 
The variable U2 discussed in section 4.2 and 4.3 are related to the a = 2 holomorphic 
differential of the genus three curve and hence to the a = 2 Lagrangian above. 

It is also interesting to explore how the second order equations following from L a 
transform under inversion of r. Dropping the space dependence, the equation of motion 
following from (|B.1|) is 

r = t- + r- ( (4 - a)r 4 -a) . (B.2) 

1 + r 4 r(l + r 4 ) V V ; ' y J 

Under a transformation R = r — t, we get 

A E>3+2a d2 

R = — — 7+ , A A (a + 2)R A + (a-2) ) . (B.3) 
1 + R 4 R(l + R A ) v v ; v ' ' K ' 

If we ignore, for the moment, the first term on the RHS, we see that the a = equation 
in ()B.2j) maps to the a = 2 term. If we neglect the velocity dependent term on the RHS, 
the a = maps back to a = 0. 

Further consider the transformation R — ~, ^ = —^5, to get 

4 d7+2q d2 

R=— — + , A A {a + A)R i + a). (B.4) 
1 + R 4 R(l + R A ) K v ; ; v ; 

If we keep only the velocity dependent term on the RHS, we find that the case a = 
from (jB.2|) maps to the case a = of (jB.4|) . 

While these transformations include relations between the a = and a = 2 La- 
grangians, they require neglecting some terms in the equation of motion so they would 
hold in special regimes where these are valid. The transformations involving r , which we 
described in section 4, give a more direct relation between the functions solving equations 
from Lq and L2. It will be interesting to find physical brane systems which directly lead 
to L a for a^O. 

C Evaluation of the solution of the collapsing D2 con- 
figuration using Weierstrass p-functions 

Here is a warm-up which makes use of the technology employed for the solution of the 
Jacobi inversion problem, in terms of Weierstrass jp-functions. We use this to recover the 
familiar result from the 2-sphere collapse. 

The integral related to the time of collapse from an initial radius r and in dimension- 
less variables is given by the expression 



Ldr=-t. (C.l) 
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After a re-scaling r = — we have 

to ro 



v/IT7 ^ dr = -t (C.2) 



i r \/l - f 4 



and 



4 
r 



Next we will perform the substitution r 2 = x ending up with 



dr . r t 

i—===u x . (C.3) 



= f ^ x (C.4) 
7i v/4a;(x 2 - 1) 

This is a definite integral over the holomorphic differential of the elliptic curve 

y 2 = 4x 3 - 4x , (C.5) 

with branch points at x — 0, ±1 and infinity. 

It is easy to see that the former integral can be decomposed into the following two 

f°° dx _ r° dx 

which are exactly the real half-period Q of the surface and the inverse of the Weierstrass 
elliptic function g 2 , g%) with invariants g<i — 4, g% — 0. Both the real and the imag- 

inary half-period (O, fi') can be calculated by contour integration after we have defined a 
homology basis on the surface. In this case take the a-cycle to be the loop surrounding 
the points (1, oo) (or equivalently, by deformation, the points (—1, 0)) and the 6-cycle the 
loop around (0, 1) and across the two sheets. The results are 

n = f - dS and n' = i f - dS , (C.7) 



with s a real, positive integration parameter. The modulus of the torus is then simply 
given by r = = i. There exists an alternative definition for this, namely r = ij^th, 
where K{k) is the complete elliptic integral of the first kind and k, k' the elliptic modulus 
and the complementary modulus respectively, with k 2 + k' 2 = 1. This implies K(k) = 
K(k') and k — k' — 1/ \[2 for the functions and integrals defined on this surface with this 
particular choice of cycles and cuts. 

From the definition of p(u) one can derive several relations between the latter and the 
Jacobi elliptic functions [221- Take a general curve 

-i, s_ f 00 dz f°° dz 

p (w) = v= = = = / = = = = (C.8) 



^/4z 3 - g 2 z - g 3 J w a/4(z - ei)(z - e 2 )(z - e 3 ) 
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and by setting 



2 9 e 2 — e 3 7 2 

7 = ei - e 3 , k = and z = e 3 + — , for e x > e 2 > e 3 , (C.9) 

e x - e 3 s 2 



1 ^ da 



the integral becomes 

v=~/ ■ (CIO) 

7 Jo v/(l-s 2 )(l-fc 2 s 2 ) 

This is simply the definition for the Jacobi elliptic integral of the first kind and is invertible 
with Sn^iW) = jv. If we use this fact and ()C.9|) we get 

P(v) = e 3 + 2 f (C.ll) 

and for the example at hand, t\ — 1, &i — 0, e 3 = —1 and the relationship becomes 

*>M = - 1+ C2 ( l 1 \ ' ( C - 12 ) 

Returning to (|C.6|) we have 

«i = fi-p -1 (a/;4,0) (C.13) 

and 



x' = p(ux-Q;4, 0) 
2 



where we have made use of the identity 

P (,±n ) = ei + (ei -;f e '- e3) , (C.15) 

- ei 

Now from ([("TH^ and by using the following properties of elliptic functions 

Cn(v,k) = - } — and Sn 2 (v, k) + Cn 2 (v, k) = 1 . (C.16) 
Cn(iv, k') 



we obtain 



x 



Cn 2 (y/2iu u -^=J 



Converting back to the original quantities x' = r' 2 = T l2 jr\ and by substituting u\ = 
—ir t/ a/1 + Tq, we recover the desired result 
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